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Abstract. Why and how may some active regions (ARs) fre- 
quently produce coronal mass ejections (CMEs)? It is one of the 
key questions to deepen our understanding of the mechanisms 
and processes of energy accumulation and sudden release in ARs 
and to improve our capability of space weather prediction. Al- 
though some case studies have been made, the question is still 
far from fully answered. This issue is now being tried to address 
statistically through an investigation of waiting times of quasi- 
homologous CMEs from super ARs in solar cycle 23. It is found 
that the waiting times of quasi-homologous CMEs have a two- 
component distribution with a separation at about 18 hours. The 
first component is a Gaussian-like distribution with a peak at 
about 7 hours, which indicates a tight physical connection be- 
tween these quasi-homologous CMEs. The likelihood of occur- 
rences of two or more CMEs faster than 1200 km s _1 from the 
same AR within 18 hours is about 20%. Furthermore, the cor- 
relation analysis among CME waiting times, CME speeds and 
CME occurrence rates reveals that these quantities arc indepen- 
dent to each other, suggesting that the perturbation by preceding 
CMEs rather than free energy input be the direct cause of quasi- 
homologous CMEs. The peak waiting time of 7 hours probably 
characterize the time scale of the growth of instabilities triggered 
by preceding CMEs. This study uncovers more clues from a sta- 
tistical perspective for us to understand quasi-homologous CMEs 
as well as CME-rich ARs. 



1 Introduction 

Magnetic free energy is thought to be the energy source of 
coronal mass ejections (CMEs). Active regions (ARs) carry 
a huge amount of free energy and therefore are the most 
probable place where CMEs come out. Lots of efforts have 
been devoted to the triggering mechanisms of CMEs. Flux 
emergence, shear motion and mass loss all could be the ini- 
tial cause of a n isol a ted CME [e.g.. \Forbes and Pnesd. ll99E ; 



Amari et all Il96rj : \Chen and ShibataV |2o6oT 



Manchester 



200.ll. No matter which one takes effect, the determina 



tive factor of the CME's launch is the force balance between 
the inner core fie l d and the outer overlying arcade s [e.g., 
I Wang and Zhand . 120071 ; \Lhl 120071 ; \SchrijveA 120091 ] . Free 
energy stored in the source reg i on wi ll be consumed when a 
CME launches [e.g.. \Sun et 



.cgio n wi ll 
aU l2012| . 



The picture of isolated CMEs is somewhat clear. How- 
ever, it is still a question how CMEs could lift successively 
in a limited region within a relatively short interval. Usually 
the energy accumulat ion is a gradua l proc e ss in time scal e 
of hours to days [e.g., \haBonte eFaH 120071 : \Li et qZ.L [ioiol ] . 
while a CME is a sudden process releasing accumulated en- 
ergy in minutes. Why and how could some ARs frequently 
produce CMEs? Does the occurrences of successive CMEs 
from the same AR mean that the source AR accumulate 
free energy quickly? The waiting time distribution of quasi- 
homologous CMEs contains clues. 

Hom ologous CMEs were defined by I Zhang and Wang] 
2002^ after t he d efinition of homologous flares 
Woodgate et all Il984| . Strictly speaking, homologous 



CMEs must originate from the same region, have similar 
morphology, and be associated with homologous flares and 
EUV dimmings. Here, we use the term 'quasi-homologous' 
to refer to successive CMEs originating from the same ARs 
within a short interval, but may have different morphology 
and associates. 

A previous study on 15 CME-rich ARs during the as- 
cending phase of the last solar cycle from 1998 to 1999 
have suggested that quasi-homologous CMEs occurred at 
a pace of about 8 hou rs, and ther e was a t most one fast 
CME within 15 hours [Chen et all |2011b| . These results 
are important for space weather prediction, and did imply 
that the accumulation rate of free energy in an AR may 
not support such frequently occurrences of quasi-homologous 
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CMEs, and the triggering mechanisms of the first and the 
following CMEs are probably different. Three scenarios were 
proposed to interpret the averagely 8-hour waiting time of 
quasi-homologous CMEs. 

Before deepening our understanding of such a phe- 
nomenon, we need to check if a similar waiting time dis- 
tribution of quasi- homologous CMEs could be obtained for 
the whole solar cycle. In this paper, we extend the period of 
interest to the whole solar cycle 23 from 1996 to 2006. In- 
stead of searching all ARs and the associated CMEs, which 
are too many to be identified manually, we investigate super 
ARs that were reported in literatures. Super ARs are those 
with larger area, stronger magnetic field and more complex 
pattern, and thought to be the representative of CME pro- 
ducers. In the following section, we present the selected data 
and the method. In Sec|3] an analysis of waiting times of 
quasi-homologous CMEs from these super ARs during the 
last solar cycle is performed. Finally, conclusions and dis- 
cussion is given in the last section. 

2 Data Preparation 

2.1 Super ARs and Associated CMEs 

Super ARs were studied by many researchers \BaA 
198lll988l:lrian et all 120021 : \Romano and ZuccarellA 1 20071 : 



numebr of super ARs almost linearly decreases with increas- 
ing CME number though there is a sharp decrease below the 
CME productivity of 3. 



Chen et all 1201 la| . It was first defined bv \Bai 1 19871 . 



as a region producing four and more major flares. In most 
studies, super ARs were selected based on several parame- 
ters, such as the largest area of sunspot group, the soft X-ray 
flare index, the 10.7 cm radio peak flux, the short-term to- 
tal solar irradiance decrease, the peak energetic proton flux, 
the geomagnetic Ap index, etc. No matter which one or 
more criteria are used, most selected super ARs are CME- 
productive (that could be seen at the last paragraph of this 
sub-section). 

In our study, we focus on super ARs during solar cycle 23. 
Instead of identifying super ARs by ourselves, we simple use 
existent lists of super ARs in literatures. To our knowledge, 
there are three lists regard ing to super ARs in solar cycle 
23. The first one is given bv \Tian et~al\ [20021 ] , who found 16 
super ARs from 1997 to 200 1 base on their selection criteria . 
The second one is given by \Romano and Zuccarelld |2007|] . 
which contains 2 6 super ARs from 2000 to 2006. The last 
one is m paper bv lCfeera et ai1|2011al ], in which 1 2 super ARs 
were id entified during the last solar cycle. Since \Chen et ali 
[2011a| | used stricter criteria, the last list is actually a subset 
of the other two. Totally, we have 37 super ARs from 1996 
to 2006. 

To identify the CMEs originating from these su- 
per ARs, we examine imaging data fro m Large Angle 
and Spectrometric Coronagraph (LASCO, \Brueckner et ali 
19951) and Extreme Ultra violet Imaging Telescope (EIT, 



Delaboudiniere et al\ Il995l ) onbard Solar and Heliospheric 



Observatory (SOHO). The identification pro cess is the 
same a s that applied bv \Wana etai\ [20 111 ) and \Chen et al\ 
[201 lb| . The CME s listed in the CDAW LASCO CME 
catalog (refe r to |http: / / cdaw.gsfc.nasa.gov / CME_list/| 
I Yashin et al\ l2004h are our candidates. Through a care- 
ful identification, it is found that a total of 285 CMEs are 
associated with these super ARs. Figure \T\ shows the distri- 
bution of the CME productivity of super ARs, in which the 
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Figure 1: Distribution of CME productivities of super ARs. 

It should be mentioned that there are 7 super ARs with 
too many large data gaps in LASCO and/or EIT observa- 
tions, and therefore their CME productivity cannot be ob- 
tained. Except them, there were 28 super ARs producing 
3 or more CMEs (called CME-rich ARs), among which 14 
super ARs generated at least 10 CMEs. The other 2 super 
ARs produced only one or two CMEs though sporadic data 
gaps existed. This fact suggests that not all of super ARs are 
CME productive. But it is de finite that super AR s are more 
likely to be CME productive. \Chen et al] |2011bf | identified 
108 ARs during 1997-1998 and found that only 14% of these 
ARs produced 3 or more CMEs. This percentage is much 
lower than that for super ARs, which is about 93% (28/30). 
In this study we focus on the 28 CME-rich ARs, which pro- 
duced 281 CMEs in total. A list of all the CMEs associ- 
ated with these CME-rich super ARs can be retrieved from 



http://space.ustc.edu.en/dreams/quasi-homologous-cmes/l 



2.2 Waiting Times 

As long as there is no large data gap, we tentatively be- 
lieve that all CMEs originating from a super AR of interest 
are recognized based on combined observations from SOHO 
LASCO and EIT. The waiting time of each CME is obtained 
according to the times of first appearance of the CME and 
its preceding one from the same super AR in the field of 
view of LASCO/C2. However, data gaps exist, and some 
CMEs may missed. If there was a large data gap between 
two CMEs from the same super ARs, the waiting time of the 
second CME cannot be obtained. Here, all data gaps less 
than 3 hours are ignored, because it is almost impossible for 
a CME to stealthily escape the field of view of LASCO in 3 
hours. 

Before analyze the waiting times of these CMEs from 
the super ARs, it has to be noted that there ar e probably 
abou t 32% of fron t side C MEs missed by SOHO \Ma et ali 
l2010l : I Wang et all l201l| . Of course, these missed CMEs 
mi ght be gene r ally we ak and faint. The statistical study 
bv I Chen et al\ [2011bl ] have suggested that the properties 
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of ARs have effects on the CME productivity, but do little 
with the kinetic properties of CMEs. Thus, it is possible 
that some CMEs originating from the super ARs are missed 
in our study, though such CMEs might be very weak and 
erupt in a gradual manner. So far. it is hard to evaluate 
how significant an influence this error will cause, and one 
may bear it in mind that the following analysis is performed 
with a bias of normal to strong CMEs. 

3 Results 

3.1 Waiting Time Distribution 

The average value of the waiting times is about f 7.8 hours. 
The waiting time distribution is shown in Figu re [21 Simi- 
lar to t hat shown in Figure 10 of the paper by \Chen et al\ 
[2011bl |. the distribution consists of two components. One 
component locates less than 18 hours and looks like a Gaus- 
sian distribution, and the other beyond 18 hours. For the 
first compon ent distribu t ion, th e peak waiting time is about 
7 hours. ln \Chen et all [2011bl |. the separation of the two 
components of the distribution is near 15 hours, and the 
first component distribution peaked near 8 hours, which are 
both slightly different than those obtained here. These slight 
differences might be caus ed by the solar cyc le variation. 

An interesting result in I Chen et al\ [2011bl ] is that any AR 
cannot produce two or more CMEs faster than 800 km s _1 
within 15 hours. In other words, the time intervals between 
fast CMEs are longer than 15 hours. If this result obtained 
during the last solar minimum also holds for the whole solar 
cycle, we could expect that any AR cannot produce two 
or more CMEs faster than a certain speed threshold within 
18 hours. However, such a speed threshold cannot be found. 
The blue line in Figure [5] shows the waiting time distribution 
for CMEs faster than 1200 km s _1 . Note that all the slower 
CMEs are ignored when we calculate waiting times for CMEs 
faster a certain speed threshold. Some fast CMEs did occur 
in the same ARs within 18 hours. For example, there were 
four CMEs from the super AR 10720 on 2005 January 15 at 
06:30 UT, 23:06 UT, on January 17 at 09:30 UT and 09:54 
UT, respectively, which were all faster than 2000 km s _1 . 
The first two CMEs were separated by about 16.6 hours, 
and the other two by about only 24 minutes. These fast 
CMEs caused ground -level enhancement (CLE) event [e.g., 
\Grechnev eFoU |2008| |. 

Although a similar result cannot be obtained, we find that 
the likelihood for an AR producing two or more fast CMEs 
within 18 hours is much smaller than normal. For all CMEs, 
68% of the waiting times are shorter than 18 hours, while 
for CMEs faster than 1200 km s _1 , the fraction decreases to 
only about 18%. The dependence of the likelihood on the 
CME speed threshold is given in Figure [3] Generally, the 
likelihood monotonically decreases as the speed threshold 
increases. When the threshold reaches to about 1200 km 
s _1 , the likelihood stops decreasing and stays between 15%- 
25%, suggesting a limit likelihood of approximate 1/5. 

The waiting time distribution for all CME s from 1999 
Febru ary to 2001 December was investigated bv \Moon et al\ 
[20031 ] . which is significantly different from the distribution 
for quasi-homologous CMEs obtained here (see Figure 1 in 
their paper). This difference rev eals that the occurrence of 
CMEs follows a Poisson process \ScaraleV Il998l : I Wheatland 
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Figure 2: Waiting time distributions for all quasi- 
homologous CMEs (black line) and for quasi-homologous 
CMEs faster than 1200 km s" 1 (blue line). 
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Figure 3: Dependence of likelihood of quasi-homologous 
CMEs occurring within 18 hours on CME speed. 



l200rj |. but that of quasi-homologous CMEs does not. In a 
statistical view, we may conclude that there are tight phys- 
ical connections between quasi-homologous CMEs, but for 
CMEs from different source regions, the connection is quite 
loose. 

3.2 Role of Free Energy Input in Causing 
Quasi-Homologous CMEs 

Sufficient free energy is a necessary condition for an AR to 
produce CMEs. Generally, the accumulation rate of free en- 
ergy could be approximately represented by the magnetic he- 
licity injection rate, which is another important parameter in 
evaluating the productivity of ARs. Magnetic helicity mea- 
sures the twists, kinks and inter- linkages of magnetic field 
lines, which indicate the complexity and non-potentiality of 
a magnetic system. The close relationship between the free 
energy and magnetic helici ty could be seen from their for- 
mulae \Kusano et al\ . 12002? ] . Thus it is not surprising that 
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a higher injection rate of magnetic helicity often implies a 
higher probability of an eruptive activity, as suggested by 
many studies [e.s. AZhana et all\200&\LaBonte et al.l. 120071 ] . 

However, it is still questionable if free energy input is a 
direct cause of quasi-homologous CMEs. Some studies did 
show that CMEs do not always occur during a quick in- 
jection of magnetic helicity or free energy, even if the stored 
free energy in an AR was much higher than that required for 
a CM E [e.g. , \Demoulin and PariaA |2009| ; I Vemareddv et all . 
2012], This issue is addressed here in a statistical perspec- 
tive from two aspects. First, we investigate the correlation 
between the CME speeds and waiting times. If free energy 
input is a direct cause, it is expected that there is some reg- 
ulation between CMEs' speeds and their waiting times, as a 
long waiting time may lead to more free energy in an AR. 
This expectation is established under the assumption that 
the injection rate of free energy or magnetic helicity varies in 
a relatively small range for different ARs, This assumption 
is statistically true bas ed on previous stu dies. For example, 
the statistical study bv \Park et al\ [20ld | suggested that the 
magnetic helicity fluxes in 378 ARs observed by SOHO/MDI 
were on the order of about 10 40 Mx h _1 , especially for those 
ARs with large magnetic flux (see, e.g., Fig.l, 3 and 4 in their 
paper). The value does not change much even if deriving 
from higher-resolution data from SDO/H MI, e.g., the helic- 
ity in jection rate in AR 11158 and 11166 \ Vemareddv et all 
l2012f| . 

Figure [4^i shows the dependence of CME speed on the 
waiting time. Overall, no clear correlation could be found 
between them, except that there is seemingly an upper limit 
in CME speed depending on the CME waiting time. How- 
ever, although the distribution is statistically true, it does 
not imply that an AR is difficult to produce a fast CME if it 
had waited too long. It is a result simply from a combination 
of two Gaussian-like distributions. The CME waiting time 
is a Gaussian-like distribution, at least for the first compo- 
nent (as shown in Figure [2| . The CME speed is actually 
also a Gaussian-like distribution. If the two quantities are 
independent, the 2-D distribution composed by them is like 
that shown in Figure [4^,. As a test, Figure [4Jd shows the 
distribution, in which the CME speeds in our sample are 
randomly associated with the CME waiting times. The two 
distributions given in Figure 2£i and are quite similar. 
It reflects that the CME speed is independent on the CME 
waiting time. 

Second, we check if the waiting time of a CME depends 
on the CME occurrence rate in the past 18 hours before its 
preceding CME. Figure[S£i shows the scattering plot between 
them. Apparently, a low or high CME occurrence rate may 
lead to a short waiting time of the next CME, and a long 
waiting time tends to appear when the CME occurrence rate 
is around 0.1 per hour. However, similar to the previous one, 
this distribution is also just a manifestation of probability, 
and contains less physical meaning. If we randomly asso- 
ciate the CME waiting times with the occurrence rates, a 
possible distribution of the data points is like that shown in 
Figure [SJd, which is statistically same as that in Figure [5^,. 
Thus the CME waiting time is independent on the previous 
CME occurrence rate. Both results suggest that free en- 
ergy input is not a direct cause of quasi-homologous CMEs 
though sufficient free energy is a necessary condition for an 
AR to produce CMEs. 



3500 
3000 
2500 
2000 
1500 
1000 
500 





o 



oeo o 







20 



40 60 80 

CME Waiting Time (hour) 



100 



' 20 



3500 



Randomized 




o o 



o 
o 



20 



40 60 80 

CME Waiting Time (hour) 



100 



" 20 



Figure 4: Upper Panel: Scattering plot of CME speeds ver- 
sus CME waiting times. Lower Panel: Same as Upper Panel, 
but the association between them is randomized. 



4 Summary and discussion 

In summary, by investigating 281 quasi-homologous CMEs 
originating from 28 CME-rich super ARs over the last solar 
cycle, we find a two-component distribution of their waiting 
times with the separation of the two components at about 
18 hours and the peak waiting time of the first component 
at about 7 hours. These results suggest a close physical con- 
nection between quasi-homologous CMEs which fall in the 
first component. Furthermore, the likelihood of occurrences 
of two or more fast CMEs within 18 hours decreases as CME 
speed increases. A limit likelihood of about 20% is reached 
when CME speed is larger than 1200 km s . 

The correlation analysis among CME waiting times, CME 
speeds and previous CME occurrence rates shows us the sta- 
tistical evidence that the free energy input is not a direct 
cause of quasi-homologous CMEs. It is well known that 
that the free energy stored in ARs may be much higher than 
that could be consumed by one single CME [e.g., \S un et all . 
2012]. Thus the direct cause of quasi-homologous CMEs is 
not the quick re-fill of free energy after preceding CMEs, but 
the perturbation by preceding CMEs, which may lower the 
threshold of eruption or trigger instabilities to cause the next 
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last two, in which quasi-homologous CMEs probably hatched 
from a long magnetic flux system or different magnetic flux 
systems in one AR. A simple/small AR should be difficult to 
frequently produce CMEs. A detailed investigation on this 
point is worthy to be carried out in future work. 
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CME. Pre-eruption flux rope is precisely balanced by out- 
ward force from i nner core field and inward force from overly- 
ing a r cades [e.g. jTorofc and K liem. 20051: \Wana and Zhantk 
120071 ; \Liil |2007( ]. A CME may reduce the constraint of its 
nearby flux rope system by removing overlying arcades, and 
cause the balance broken. A s shown in the numerical sim- 
ulation by I Tor ok et~al\ [201 ll ]. which was designed to study 
the physical mecha nism of a global sympathe tic eruptions 
on 2010 August 1 \Schrijver and TitleY 1201 1| . the second 
and third eruptions were actually caused by preceding erup- 
tions. In their eruption processes, the preceding eruption 
caused the overlying arcades reduced through reconnection, 
and then instability devel oped. A similar resul t was ob- 
tained in the simulation bv \Bemvorad et~al\ [20121 ] . in which 
the second CME was caused by the rearrangement of coronal 
magnetic field after the first CME. 

Connecting the above picture to the peak waiting time 
of 7 hours, we may speculate that the 7-hour waiting time 
probably characterizes the average tim e scale of the growth 
of instabilities. In our previous work \Chen et all 1201 lbl ] . 
we proposed three scenarios to interpret the peak waiting 
time. Here we may tentatively narrow down them to the 
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